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Abstract 

We study the dynamics of ultracold attractive atoms in a weakly linked two potential wells. We 
consider an unbalanced initial state and monitor dynamics of the population difference between the 
two wells. The average imbalance between wells undergoes damped oscillations, like in a classical 
counterpart, but then it revives almost to the initial value. We explain in details the whole behavior 
using three different models of the system. Furthermore we investigate the sensitivity of the revivals 
on the decoherence caused by one- and three-body losses. We include the dissipative processes using 
appropriate master equations and solve them using the stochastic wave approximation method. 

PACS numbers: 03.75.Gg, 03.75.Kk, 67.85.Dc, 05.30.Jp 
Keywords: 
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I. INTRODUCTION 



In 1962 Josephson predicted the current flow associated with coherent quantum tunneling 
of Cooper pairs through a barrier between two superconducting electrodes. This device is 
called a Josephson junction and the flow is proportional to the sine of the phase difference. 
Jospehson also predicted that the time derivative of the phase difference is proportional to 
the voltage across the barrier. The corresponding relations can be found in l|-|3|. The first 
superfluid Josephson junctions were realized in superfluid 3 He in 1997 j4-6]. The Josephson 
junction dynamics relies on the existence of two coupled macroscopic quantum states. 

With the advent of Bose Einstein condensates of weakly interacting gases a new experi- 
mental system has become available for the quantitative investigation of Josephson effects 
in a very well controlled environment [zj|. Note that the Josephson junction in this system 
consists of the two localized matter wave packets in the two wells coupled via tunneling of 
particles through the potential barrier. The authors of jsj presented the experimental real- 
ization of the atomic Josephson Junction and compared the data obtained experimentally 
with predictions of a many-body two-mode model and a mean-field description. 

In our paper we study two mode Bose-Hubbard model with the attractive interactions. 
We focus on the dynamics of the state that is initially in a condensate state (all atoms 
occupying the same single particle state) with nonzero population imbalance. During the 
time evolution there are Rabi oscillations in the mean value of the population imbalance. 
These oscillations exhibit collapses and revivals, analogous to the well known phenomena 

IT 

in quantum optics (see Jaynes-Cummings model [9J). The revivals observed in quantum 
optics have been viewed as a proof of the existence of photons i.e. one could not explain 
this phenomenon within classical electromagnetic field interacting with atom. In this study, 
we explain the mechanism of the revivals in the double-well system in the same manner as 
revival predicted in James-Cummings model. However, we do not claim that the revivals 
are the evidence of particle existence. On the contrary, we reconstruct the revivals using 
only concepts from classical physics. Our classical model also predicts correctly collapses 
and Rabi oscillations. Next we include in our consideration particle losses. Our numerical 
calculations prove that the revivals are extremely sensitive to decoherence - loosing just a 
few particles destroys it almost completely. 

The paper is organized as follows. In section [TT] we introduce a two mode Bose Hub- 
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bard model. We focus on the attractive interaction case and shortly derive the continuous 
approximation. This approximation maps the two mode Hamiltonian onto the Hamilto- 
nian of a fictitious particle. We obtain a one- dimensional Schroedinger equation with an 
effective Planck constant h e ff inversely proportional to number of particles. The shape of 
the potential determines the dynamics of the system. This approach offers a better insight 
into the dynamics of the system. In section II III we explain the origin of Rabi oscillations, 
their collapses and revivals. Furthermore we compute the timescales characteristic for these 
phenomena. The section IIVI is devoted to a decoherence in the system and we discuss 
its influence on the dynamics. We include one- and three-body losses in the system and 
investigate their role in damping of revivals. 

II. MODELS 

We consider an evolution of the gaseous attractive Bose-Einstein condensate in a double- 
well potential in the two mode approximation. In the limit of relatively weak interaction 
the system can be described by the Bose-Hubbard Hamiltonian 



where J is the tunneling energy, U is an average interaction energy per pair of atoms, N is 
a total number of particles and operators a\ and a\ create one atom in left and right well 
respectively. In all cases except these with particle losses, it is useful to omit constant terms 
in ([!]) and rewrite it in the dimensionless form 



where 7 = ^j-. This form of Bose-Hubbard Hamiltonian determines our dimensionless time 
unit r :— M. To get a better insight in the structure of the Hamiltonian (more 'intuitive 
one') one can use a continuous approximation, valid for large N. Consider the expansion of 
the state of system in the Fock basis: \ip) = Yln=o' l Pn\ n ' N — n), where the ket \n, N — n) 
means a state with n atoms in the 'right' well and rest of them in the 'left' well. The 




(1) 
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Schrodinger equations for the coefficients ip n can be written as 

= -Jn ( 2n ~ N ?^n - y/n{N-n + l)^n-i 

-V(n+1) (iV-n)^ n+1 . (3) 

Assuming that iV ^> 1 and ip n are negligible near the 'boundary' (for n close to zero and N) 
we can approximate equation (j3J) getting 



where z n = — n . In the same limit a discrete set of the coefficients i/j n may be replaced 
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by a 'wave function' ip(z). Upon replacing (ij) n -i + ip n +i — 20 n ) — jv* JjaVK^) the equation 
(HI) takes the form 



where 



l Nfr= —d^ + V{z) ^ (5) 



V(z) = -y/T^- |z 2 . (6) 



We further simplify the equation fl5]) through omitting the term a/1 — z 2 in front of the 
second derivative: 

The latter equation describes a fictitious particle of mass m e ff = | with the effective Planck 
constant h e ff = j^. We will consider only evolution limited to small mean value of position 
(z). In these cases replacing the effective mass A ^^_ z -2 from Eqn. [5] with only m e // = \ 
will not change qualitatively evolution. The shape of the effective potential is shown in Fig. 
[TJ With increasing 7 the shape of the effective potential ([6]) turns from a single minimum 
located at z = 0, to a double- well structure for 7 > 1. 

Let us define the family of initial states, which we will consider in the rest of the paper. 
As an initial state we take a coherent state with iV-atoms in the same single-atom state 

m) = I + |0>. (8) 



zo = <0(O)|i|V(O)>, (9) 



f Nl 

The zq is the initial imbalance 
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FIG. 1: (Color online) Shapes of effective potential for various values of 7. 
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FIG. 2: Numerical solution of Bose-Hubbard model. Average relative atoms' imbalance between 
wells (z) for 1000 particles, the initial average imbalance zq = 0.4, the initial phase 4>q = and 
7 = 1.0. 



where 



z := (a[ai — a/ 2 a2)/N 



(10) 



is a population imbalance operator. The parameter O can be interpreted as a relative phase 
between two wells. 

In the continuous approximation the initial state corresponding to the state (IE]) can be 
constructed upon it's expansion in the Fock basis 



] j^p n/2 q {N - n)/2 e i{N - n) ^\n, N-n), 



(11) 



where p = (^Sp-) and q — 1 — p. The square of modulus of the expansion coefficients in 
(ITT]) is just a probability function for the binomial distribution. In the continuous limit ( 
N — > oo and z n — > z) they converge to a Gaussian function and we get 



The above discussion is presented in more details in [lOj, 111]. 

The continuous model together with the initial state proposed above give us a clear 
picture of what will happen in the system. We have a Gaussian wave packet moving in the 
potential V(z). The shape of the potential will be crucial for the dynamics. In all further 
considerations we focus on the case 7 = 1, where in the effective potential the harmonic 
terms z 2 are canceling and we approximate it by the first two leading terms: 

V ef = / « "I + J- (13) 

Having chosen 7 and the initial state |^(0)) the only free parameter is the number of particles 
N. 

At last, let us focus on the dynamics. We evolve the state \4>(t)) using the Schroedinger 
equation ([3]). Fig. [2] shows the (ip(t)\z\ip(t)) in the case iV = 500. Here we observe Rabi 
oscillations, with distinct collapse at time r ~ 150 and the revival at r ~ 750. 

In the next sections we focus our attention on the variable (%lj{t)\z\^{t)) . 

III. THE ORIGIN OF RABI OSCILLATIONS, COLLAPSES AND REVIVALS. 
A. Rabi oscillations 

Due to the tunneling between the wells for most initial conditions (z) will oscillate peri- 
odically in time. This phenomena can be easily illustrated within the classical field theory 
(substituting the annihilation and creation operators by complex numbers). This approxi- 
mation may be understood as a transition from the Schrodinger equation ([7]) to the Newton 
equation 

d 2 z dV 

= -m (14) 

The motion of the fictitious particle described by the above equation is periodic. In the 
spirit of quantum optics we call these the Rabi oscillations. Using the Newton equation we 




t [ T Ra bl ] 



FIG. 3: (Color online) Comparison between evolution of (z) given by Bose-Hubbard model for 
2000 particles (red - dashed line), for 8000 particles (green, dot-dashed line) and classical model 
(black, solid line) for zq = 0.4, </>o = and 7 = 1. The scale on x-axis corresponds to analytical 
calculated time of the Rabi oscillation. 



calculate their period Tj^. The formula following from Eq. (1141) reads 

dz 



■ Rabi 



y/2(E-V(z))/m eff 



(15) 



The energy E is just a sum of a potential and a kinetic energy ^fr^C?- From the initial 



state it's clear that E = V(zq) 



2 dt' 2 ■ 

z 4 

y. Substituting this value into f fT5"j) we can estimate 

_ 4y^r (I) 

In Fig. [3] we compare this result with numerical simulation of Bose-Hubbard model for a 
short time evolution getting a good qualitative agreement. From this figure we can also 
conclude, that the approximation become accurate in the limit N — > 00. 



B. Collapse 

The single trajectory of Newton equation contains oscillations which continue forever, 
so it does not describe the phenomenon of a collapse. We argue that one can obtain the 
collapse using Newton equation with a set of random initial conditions, with appropriate 
jrobability distribution. Our method will be similar to a semiclassical model exploited in 



12l | . This approach is also equivalent to a stochastic classical field theory 
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FIG. 4: (Color online) Comparison between evolution of (z) given by Bose-Hubbard model (black 
- solid line) and semiclassical model (red - dashed line). 

Central issue is: is there an appropriate probability distribution corresponding to the 
initial state f lT2|) . One of the candidates is the Wigner distribution if it is positive. Since the 
classical equation we use is the Newton equation corresponding to the Schrodinger equation 
([7]) the Wigner function can be calculated as 



The procedure used to describe the dynamics will consist of the following steps. We draw a 
set of (in our case 50000) initial conditions, propagate them for some interval of time t and 
form the distribution again. Using this new distribution we can calculate any expectation 
value. In particular in Fig. [4] we present the evolution of (z) obtained this way and compare 
it with the numerical solution of the quantum Bose Hubbard model. The agreement is 
satisfactory and gives a correct prediction of the collapse. It indicates that the collapse is 
a result of initial conditions uncertainty (spread). We clearly see that the revivals are not 
present in this approach. This is consistent with the understanding of this phenomena in 
quantum optics. 

In Fig. [5] we present the reconstruction of the classical distribution in the phase space 
z — cj) at different stages of the evolution. This picture illustrates generation of the damping 
of Rabi oscillations. When the distribution reaches a shape of ring, the average (z) becomes 





(17) 
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(a) r = 



(b) r = 40 




FIG. 5: (Color online) The evolution of classical distribution in z — <j> space, initially equal to the 
Wigner function. The bright regions corresponds to high probability of finding an atom in this 
region. Parameters in simulations N = 2000, zO = 0.4, <j)Q = 0.0 and 7 = 1. 

and the collapse occurs. The merging of the distribution's tails is due to the dispersion of 
energy in the initial conditions 

AE = yj (H*) - (H)* = l/VN + 0(±), (18) 

and the dependence between the time of one period and the energy Trom oc E~ x ^ (see Eq. 
(fl5]) ). Using (fl~8l) and (fl~5l) one can show that 

T co i oc VN. 

In Fig. |6]we compare our prediction with the numerical results getting a good agreement. 
The derivation presented above is only a very rough estimation of the real process. In fact the 
collapse occurs when the distribution 'winds up' three or four times so when the distribution 
smears uniformly over ellipses in the phase space (Fig. [5]). However we think that the 
presented mechanism of collapse is intuitive and useful for the qualitative predictions. 



C. Time of revival 



The most famous revival is the one in the evolution of a two level atom interacting with 
state of a single electromagnetic mode. This was first predicted in js], and finally observed 
13[. The theoretical considerations used to calculate the revival time in 9[ may be also 
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FIG. 6: Timescales of collapses for different number of atoms N - results of numerical simulation 
of Bose-Hubbard model (black circles). As predicted in the text the dependence is like y/~N (solid 
line - fitted curve). 
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FIG. 7: (Color online) Comparison between the numerical solution of Bose-Hubbard model (solid, 
black line) and Schrodinger like equation (dashed, red line) for N = 2000 particles and initial 
imbalance zq = 0.4. 

applied to our system. However, to follow the derivation we need the structure of energy 
levels. 

On the other hand we see in Fig. [7] the Schrodinger equation reveals the revivals. The 
time of revival in the frame of Schrodinger equation agree qualitatively with T rev from Bose- 
Hubbard model. Thus we compute the time approximating the system with the Schrodinger 
equation. In this way we gain the knowledge of energy levels, which for particle moving in 



11 



a smooth potential may be easily estimated from WKB theory 

We decompose the state of the system \ip{r)) in the basis of the eigenstates 

| V(t)> = Z n A n (0)e iE ^ff | E n ) , 



(19) 



where \E n ) denotes the n th eigenstate of the Hamiltonian ([7]). Thus, the evolution of (z) is 
given by the equation: 



(z)(r) = i: n , m A n (0)A* m (0)e- i ^- E -^ff(E m \ z \ E n ). 



(20) 



Here, we estimate the energy spectrum using the Wentzl-Kramers-Brillouin approxima- 
tion (WKB, see for instance 14J). The stationary version of the Schrodinger equation (J7J) 
has a form: 

d 2 A 



dz- 



+ K 2 (z)lp n = 0, 



(21) 



where k(z) - 
ing relation 



2m e ff En . i v ^ . One of conclusions of the WKB approximation is the follow- 

n e/f 



(22) 



which is widely used to estimate the eigenenergies E n . The points z m in and the 
turning points for a particle with energy E n moving in the potential V(z). The integral in 
( 122]) has an analytical solution, 

~ \3/4 



k(z) dz 



N (E n ) (1/4) 



16 T (7/4) 

where E n = 8(E n + 1). Hence we approximate the eigenenergies using the following formula: 

4 

D 



E„ 



N 



(2ra + 1) 



(23) 



Here D = 8 ^nfl)^ • ^ e define the number n, which corresponds to the index of state with 
the energy equal to the average energy in the system. It is given by 



N . 1 
" = 2^°-2 



(24) 



The investigated sum (1201 will be maximal if all the components have the same phase. 
Obviously this condition couldn't be fulfilled for all terms in ( 120|) . but it is sufficient to 
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FIG. 8: Times of revival for different number of atoms N. Dots - the numerical simulation of the 
Schrodinger equation (zq = 0.4,7 = 1), solid line - the analytical estimation. 

ensure it for major terms. So let us investigate the differences between phases for terms 
with indices close to n. Thus, we are looking for a time T rev that satisfies the relation: 



After substituting equations fl24"|) and fl23l) to the condition fl25|) we get the following esti- 
mation for the revival time 



The test of the linear dependence T rev with respect to the number of particles iV is 
presented in Fig. El 

D. Classical revival 

Here we show, that though we start from the state that entangles atoms from both poten- 
tial wells OH]) the main features (including revivals) of the quantity (z)(t) can be reconstructed 
classically. More specifically we find some classical distribution in the z — <j) space and we 
use it to distribute initial conditions fot the Newton equation ( jl"4"l) . After averaging these 
classical evolutions we get some average trajectory (z)(t). Here, we present our method in 
the system of N = 2000 atoms and initial imbalance zq = 0.4. 

From Fig. [2] we see that the revival occurs about r = 740. We will find such initial 
conditions for Newton equation, that the revival will be enforced. One way is to choose such 



Trev (E n +2 — 2E n+ i + E n ) — 2irh e ff. 



(25) 




(26) 
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FIG. 9: (Color online) Comparison between the numerical solution of Bose-Hubbard model (solid, 
black line) and average over classical trajectories(dashed, red line) for 2000 particles. 



discrete set Zq that T rev will be a multiple of the Rabi periods TR a u(zQ l> ): 



T 



n G N 



» 



A^nY (|) 



(27) 



J-RabiyZo ) J- rev 1 - \ 4 J 

Thus when we evolve the point z^ using the Netwon equation, at time T rev it has to return 
to its initial position. However, to satisfactory reconstruct the mean value (z) we have to 
compute average of classical trajectories using some weights. We propose weights based on 
the Wigner distribution: 



(n) (n+l) 
~0 +z 



win] 



W(z,(j)) d(j)dz 



Thus our final approximation is 



; (t) = X>w^ (n) w, 



n=0 



An) 



(28) 



defined in 



where z^ (t) is the solution of Newton equation with initial conditions Zq 
( !27l) and % \t=o— 0. The comparison between our approximation f j28|) and the numerical 
solution of the Bose-Hubbard model is presented in Fig. [9j Such constructed solution ( 1281) 
reveals not only the times of revivals, but also the time of collapse and the period of Rabi 
oscillation. 
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IV. DECOHERENCE 



So far most experiments with BEC in a double-well potential were performed with atoms 
with repulsive interactions js), Il6 |. On the other hand experimental data can be collected 
over very long time, up to a few seconds. Assuming that the similar experiments will be 
possible in the near future, what will be the biggest challenge on the way to observe revivals? 
We found that one of the crucial parameters is the number of atoms. Hence we think that 
it is very instructive and important to study our system in the context of the particle loss. 

In this paper we demonstrate two mechanisms of the decoherence. The first of them is 
due to collisions with a hot background gas. A typical effect of such a collision is a loss of a 
cold atom from the trap. The second process under consideration is a three-body loss due 
to a recombination. If three condensed atoms meet, two of them may create a bound state 
and push the third to fly away carrying the energy gained from bound atoms. In the next 
subsections we review some of the main features of these decoherence sources and present 
master equations which cover both the free evolution of BEC and atoms' losses. These 
master equations are of the type 



dp i 
~dt ~ h 



Lp, (29) 



P,H 

where H is the Hamiltonian of a quantum subsystem and the superoperator L has the form: 

Lp = felpCk - C\C kP - pC\C k ) . (30) 

We solve the equations numerically in the frame of the quantum trajectory approximation 

In this 



(also called the stochastic wave approximation) described in details in [17 
approximation the density operator of the system is written as the average of the stochastic 
wave functions 

1 - 

= J m m E \^o Ch (tm?to Ch (t)i (si) 

rn=l 

where \ty™ t odSf)) * s a m_ th stochastic wave function at the time t. The function \ip^ och (t)) is 
also called a single quantum trajectory. It is obtained in the stochastic evolution described 
in the following. 

In the evolution of the system during some small interval At one of the elementary 
processes occurs: either atoms undergo the counterpart of free evolution or one of the 
quantum jumps. 
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The quantum jumps are just actions of the operators C k on the state of the system 



\^ stock) H> C k \lp stock)- 



(32) 



The k-th quantum jumps occurs with the probability pt = (ClCk)At, where the average 
(C k Cf.) is computed for the stochastic wave and generally it differs between trajectories. 
The 'free evolution' is represented by the effective operator 



H = H - ^C\C k . 



(33) 



If none of the quantum jumps occur (what happens with the probability 1 — YlkPk) the 
system undergoes the transformation 



\ipstoch) !-> e '» At \i) s toch)- 



(34) 



The operators H and C k are non-hermitian, so the stochastic wave function must be 
normalized after every time step At. 

A single stochastic trajectory is sometimes interpreted as a single experimental realization 
of the quantum evolution. This was the case of the famous experiments showing the quantum 



jumps 
in BEC 



20 



he stochastic evolution was successfully app 



211 ] and the interference of two Fock states 22J. Recently, the quantum jumps 



ied in the study of the decoherence 



trajectories (generated in the frame of the large deviation method) turned out to be very 



useful for studying the phase transitions 



231. 



A. One-body losses 



Due to non-perfect vacuum all experiments with ultracold atoms are performed in the 
presence of a background residual gas. Atoms from the residual gas affect the condensate via 
collisions with ultracold atoms. We assume that the typical energy of an incoming atom from 
the background is so high that during the collision the condensed atom is kicked out from 
the double-well trap. This kind of losses plays the major role for long lasting experiments, 
which may be the case of presented revivals. We propose the following phenomenological 
master equation which includes the one-body losses 

2 



dtp 



h 



H,p 



+ Kx 



k=X 



+ 



(35) 
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FIG. 10: Evolution of (z)(t) (a) without losses, (b) with in average one quantum jump before 
revival, (c) with in average four quantum jumps before revival, (d) Single realization of stochastic 
wave approximation. The phase of Fock components of stochastic wave change randomly. This 
process leads to jumps of average (z), which significantly differs from its average evolution. In the 
all figures zq = 0.4, the initial total number of atoms N = 1000 and 7 = 1. 

where H is given by equation (CQ) and k\ is a rate of one-body losses. In a case of losses 
caused by the residual gas, the parameter k± is just an inverse of the lifetime in a trap 24]. 
The form of equation (135)) has been proven in the case of suppressed tunneling (J = 0) in 



Ref. 



20]. 



In the case of the master equation (135]) the evolution operator and quantum jump oper- 
ators have the form 

H = H 

C 2 = V2^~ia 2 . (36) 

Note that the free evolution operator H is given just by the hermitian operator H. According 
to the equation ( 133]) this operator should have another form, namely 



Hi = H — ihKi(a[di + = H — ihK\N, 

where N is the operator of the total number of atoms in both wells. Then, after the 
nonhermitian evolution 

Wstoch)(i~ + Ar) = e » T \^ s toch)(r), 
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the stochastic wave function \^ s toch){j + At) ought to be normalized. However, in this case 
the normalization would just mean division the state by exp(— Nk\At), where N would 
be the total temporary number of atoms. Thus the procedure of evolution with the non- 
hermitian operator and then the state's normalization is equivalent to the free hermitian 
evolution, resulting from operator H defined in Eqn. ([I]). 

We carried out the series of stochastic evolutions for different damping parameters k±. 
In Fig. [10] we show the evolutions with K\ — (Fig. flDT a). with such damping that in the 
system up to the first revival one particle is lost (Fig. fTDTb) and four particles are lost (Fig. 
[TUb). We obtained these results averaging over only 100 quantum trajectories. Figures ITOT d) 
displays just a single quantum trajectory. We see a fast damping of revivals, even if the loss 
of the total number of atoms is relatively small. Even for just four lost atoms, the amplitude 
of revival decreases more than twice. It gives us the following qualitative criterion for the 
upper value of parameter K\, 

n\Nt rev < 1. 

Assuming number of atoms N ~ 1000 and t rev ~ 50fis (estimated from typical parameters of 
double- well potentials) we get K\ < 20/s. In typical experiments with BEC this parameter 
is two orders of magnitude smaller 24[]. Thus, the mechanism of one-body losses shouldn't 
prevent the revivals. On the other hand we haven't included other sources of one-body 
losses, like for instance due to heating. 

One may observe, that the quantum trajectories (Fig. [TOT d) differ significantly from 
the result of averaging. Even though in a single stochastic realization the amplitudes of 
revival are typically as large as in the case without losses. However, the decrease of revivals' 
amplitude in average is fairly intuitive. Let us remind, that the time of revival is proportional 
to the number of atoms. If we permit atomic losses the number of particles may change 
so the time of revival occurrence become shorter. The average (z) has solution oscillating 
in time with period Tr^. Thus the averaging over many stochastic trajectories is just a 
destructive interference, which results in the revival's damping. 



B. Three-body losses 



The main source of three-b ody losses is the recombination of condensed atoms, discussed 
in details in both theoretical 25 



26| and experimental (271 . |28| works. Our study of three- 
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FIG. 11: Evolution of (z)(t) (a) without losses, (b) with in average one quantum jump before 
revival, (c) with in average four quantum jumps before revival, (d) Single realization of stochastic 
wave approximation. In the all figures zq = 0.4, <j)Q = 0, the initial total number of atoms N = 1000 
and 7 = 1. 

body losses is just a simple extension of the previous subsection. We propose, also in 
phenomenological manner, the following master equation 

% r / A3 

h 



dp 
~dt 



Ha,P 



k=l 



{a k ) 3 ,p(al 



(a fc ) 3 p, (d\ 



(37) 



where k% is in the relation to the recombination event rate constant K 3 (see for instance 



20 



29 



3Cj) 



^3 



K3 

6 



»| 6 d 3 r. 



Like in the previous case - the equation above was derived in the limit of suppressed tunneling 
between wells. The results of the master equation in our case, namely J > 0, we treat only 
as a qualitative estimation. 

The non-hermitian evolution operator and quantum jumps operator have forms 



H 

c 2 



H — ihn? 



v / 2K 3 "a 2 . 



(38) 



We present results of numerical simulations in Fig. [TTJ Figures (a), (b) and (c) display 
the evolution up to the first revival without losses, with on average one quantum jump 
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and with on average four quantum jumps, respectively. Like in the previous subsection the 
revivals are strongly suppressed even with relatively small losses and a single realization 
(Fig. ITTT d) of stochastic process differ from the average evolution. 

Assuming that to observe the phenomenon, losses should be smaller than one atom until 
the revival, we get the qualitative criterion for k 3 



K 3 NH rev < 1. 



For example for N = 1000 and t rev = 50/xs one get k 3 < 0.5 x 1C 



-10 



/s. From the rough 



estimation of the possible experimental setups (K3 < 10 _28 cm 3 /s ( 30|, |3l|) and the size of 
each BEC in the wells about few micrometers) we compute the parameter K3 of the same 
order of magnitude. 



V. CONCLUSIONS 



In this paper we investigated attractive BEC in a double-well potential. We first recog- 
nized three essential time scales characteristic for the evolution of the relative number of 
atoms' imbalance between the two potential wells. The first timescale, the Rabi oscillation, 
has a purely classical origin. It can be understand within the frame of the mean field ap- 
proximation. The second timescale is the time of collapse. We have explained it within the 
semiclassical picture, using differential equation from the mean field approximation but for 
non-deterministic initial values. Eventually we explain the phenomenon of revivals - this 
turns out to be exactly the same as for revival in the Jaynes-Cummings model. However 
we also show that the phenomenon does not require the quantum model. It may be also 
understood using theoretical tools of classical physics. This means that these revivals are 
not the evidence of quantum correlation in the system. 

We investigated also sensitivity of revivals for one- and three-body losses. We showed 
that even if the fraction of the lost number is small (i. e. 0.1%), the revival is strongly 
suppressed. On the other hand the rate of one-body losses in the recent experiments is so 
small, that this kind of losses should be negligible. On the contrary, the three-body losses 
might be an important limitation. We do not exclude the revivals from being measured but 
our results indicate that experiments should be prepared for relatively low densities. 
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